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Abstract 



f^*) ' The results of three-dimensional MHD numerical simulations are used to investigate the 

characteristic properties of the magnetic-field structures in the accretion disks of semi- 
detached binary systems. It is assumed that the intrinsic magnetic field of the accretor star 

P^ ' is dipolar. Turbulent diffusion of the magnetic field in the disk is taken into account. The 

SS Cyg system is considered as an example. The results of the numerical simulations show 
the intense generation of a predominantly toroidal magnetic field in the accretion disk. 
Magnetic zones with well defined structures for the toroidal magnetic field form in the disk, 
which are separated by current sheets in which there is magnetic reconnection and current 

^j i 

dissipation. Possible observational manifestations of such structures are discussed. It is 
shown that the interaction of a spiral precessional wave with the accretor's magnetosphere 
could lead to quasi-periodic oscillations of the accretion rate. 

1 Introduction 

m 

There are a large number of close binary systems in which the magnetic field plays an important 
role in mass transfer and accretion. First and foremost, such systems include magnetic and 
intermediate polars, as well as X-ray binaries [UEIE]. These are semi-detached binary systems 
in which one of the components (the donor) fills its Roche lobe, allowing material from its 
envelope to overflow to the other component (the accretor) through the inner Lagrange point L\. 
The accretors in magnetic and intermediate polars are white dwarfs with surface magnetic fields 
of 10 7 -10 8 G and 10 4 -10 6 G, respectively. The accretors in X-ray binaries are neutron stars with 
intrinsic magnetic fields of 10 12 -10 13 G. It is believed that mass transfer in intermediate polars 
and X-ray binaries can lead to the formation of accretion disks around the compact objects 
[TJ [2] . The intrinsic magnetic field can appreciably influence the structure of the accretion disk 
and determine the character of accretion onto the star [U EJ El E]- 

In our previous studies [8j |9], we presented results of numerical simulations of the structure 
of MHD flows in semi-detached binary systems. These showed that the main parameters of 
the accretion disk, such as the accretion rate and characteristic density, can change when the 
magnetic field is taken into account. However, the structure of the magnetic field in the disk 
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was not considered. The magnetic field in the accretion disk around a compact object can 
be strengthened due to differential rotation, radial motions and dynamo action. Diffusion, 
turbulent dissipation, and magnetic buoyancy can weaken the magnetic field. The action of 
these effects can lead to the formation of a magnetic field with a complex structure, since 
different effects may dominate in different parts of the disk. Note that, in contrast to individual 
objects (a star, accretion disk, etc.), an important role can be played in close binaries by specific 
mechanisms for the generation of magnetic field [TO]. For example, the Herzenberg dynamo may 
operate in magnetic polars [HI H2] , leading to the generation of magnetic field in the envelope 
of the secondary. It is believed that the decelerating electromagnetic torque that arises in this 
case leads to synchronization of the rotation of the components in these systems (2] . 

The generation of toroidal magnetic field due to differential rotation dominates in the inner 
part of the disk [2j. The character of the generated field is determined by the rotation law in 
the disk. However, effects associated with the presence of a poloidal velocity in the disk, which 
redistribute the magnetic field throughout the disk, may be important here p21 HH [15] . The 
magnetic field can interact with waves arising in the inner regions of the disk [161 E] > leading 
to quasi-periodic variations in the accretion rate onto the star [9]. 

In the outer part of the accretion disk, magnetic field may partially be generated by the 
dynamo mechanism. Both a laminar dynamo (due to non-axially symmetric motions) [IB] , and 
a turbulent au dynamo [TDJ EDI EI] may operate in accretion disks in close binary systems. The 
dynamo generation of magnetic field in accretion disks has been considered, for example, in 
[22| |23| EH ESI El EE] . The dynamo generation of magnetic field requires that the mean helicity a 
for non-axially symmetric and turbulent motions of the gas not display mirror symmetry about 
the equatorial plane of the disk. Physically, this means that the numbers of right-handed and 
left-handed vortices are not equal, since the Coriolis force creates additional winding of the 
vortices. When a > the quadrupolar mode of the magnetic field is dominant above the 
equatorial plane [27] . However, detailed numerical simulations of magnetic turbulence [281 ES] 
due to the development of magnetorotational instability in accretion disks [301 EI] show that 
the opposite situation can also be realized under certain conditions [221 E3]- In this case, the 
dipolar component of the magnetic field becomes dominant. Note that this component of the 
field is preferred for the generation of bipolar outflows from the disk via the magnetocentrifugal 
mechanism [34J. 

Here, we use three-dimensional numerical simulations to study the structures of magnetic 
fields in accretion disks in semi-detached binary systems. We consider the SS Cyg system as 
an example, whose white dwarf has a dipolar magnetic field. The numerical models take into 
account radiative heating and cooling, as well as diffusion of the magnetic field due to current 
dissipation in turbulent vortices and magnetic buoyancy. 

This paper is organized as follows. Section 2 briefly describes the model and numerical 
method used. Section 3 presents the results of the numerical simulations. In Section 4, we 
consider a simple analytical model for magnetic-field generation in accretion disks in semi- 
detached binary systems. Our main results are summarized in the Conclusion. 

2 The model 

We will describe the flows in a semi-detached binary using a Cartesian coordinate system (x, 
y, z) in a frame rotating with angular velocity CI about the common center of mass. The 
coordinate origin is specified to be the center of the accretor, and the center of the donor is 



located a distance A from the accretor on the x axis. The z axis is directed along the rotational 
axis of the system. 

The MHD equations can be written 

^ + V-(pv) = 0, (1) 

dv 11 

— + (v ■ V) v = — VP --(Bx(VxB)) + 2(vxO)- V$, (2) 

f)T2 

— = Vx(vxB-r / (VxB)), (3) 

pT^ + (v-V)^=n 2 (r-A) + £-(VxB) 2 . (4) 

Here, p is the density, v the velocity, P the pressure, B the magnetic field, s the entropy per unit 
mass of gas, n = p/m p the number density, m p the proton mass, 77 the magnetic viscosity, and $ 
the Roche potential. The radiative heating and cooling T and A have complex dependences on 
the temperature T [3S1ESIEZIEH])- Our numerical model uses a linear approximation for these 
functions in the vicinity of the equilibrium temperature T = 11230 K [SHUHHO], corresponding 
to an effective temperature for the accretor of 37000 K. The term 2 (v x O) in the equation 
of motion (J2J) describes the Coriolis force. The last term on the right-hand side of the entropy 
equation @ describes heating of the matter due to current dissipation. The density, entropy, 
and pressure are related by the ideal-gas equation of state, s = cy ln(P/p 7 ), where cy is the 
specific heat of the gas at constant volume and 7 = 5/3 is the adiabatic index. 

We consider here a semi-detached binary system with a forming accretion disk and with 
synchronous rotation of the components, when the rotational period of the accretor is equal to 
the orbital period. We will also assume that the magnetic field of the accretor is dipolar, so 
that the field strength is determined by 

B = 3(m-r)r _ m 

where m is the magnetic moment of the accretor. The direction of the magnetic moment 
m need not coincide with the direction of the angular velocity of the binary system O. The 
components of m are m x = m sin 9 cos 0, m y = m sin 9 sin 0, and m z = m cos 9, where m is the 
magnitude of m, 9 the angle between m and the z axis, and <j) the angle between the x axis 
and the projection of m onto the xy plane. 

Our numerical modeling takes into account the effects of magnetic-field diffusion [HI [9]. 
Turbulent diffusion of the magnetic fields in accretion disks is determined by two main effects 
[2]. The first is magnetic reconnection and dissipation of currents in turbulent vortices. The 
second is associated with buoyancy of flux tubes of the toroidal magnetic field generated in 
the disk as a result of differential rotation. Both these effects are included in our model. Note 
that the magnetic viscosity coefficient due to magnetic buoyancy depends on the strength of 
the generated magnetic field. Therefore, the turbulent diffusion of magnetic field in the disk as 
a whole is non-linear. 

We used the Nurgush three-dimensional parallel numerical code for the numerical simula- 
tions [10] • This code is based on a higher-order, Gudonov-type, finite-difference scheme. The 
original unified-variable technique for the MHD equations [JT] makes it possible to use an 
adaptive grid in the numerical code. Our computations used a geometrically adaptive grid that 



became denser in the equatorial plane and toward the surface of the accretor. This enabled 
us to appreciably increase the resolution of the vertical structure of the accretion disk and in 
the region of the accretor magnetosphere. In order to minimize the numerical errors in the 
difference scheme, we computed only the magnetic field induced by currents in the accretion 
disk and in the outer envelope [12~| H3] . We used the eight-wave method [El HI] to clean the 
divergence of the magnetic field. 

The equations describing the diffusion of the magnetic field are non-linear. Therefore, the 
application of explicit methods to solve for this field would lead to too rigid constraints on 
the time step used. In our approach, this equation is solved numerically using an implicit, 
locally one- dimensional method with a factorized operator [45] . Another difficulty is that this 
equation contains mixed derivatives in the curvilinear, non-orthogonal coordinate system due 
to the geometrically adaptive grid used. Our method solves this problem via regularization 
of the factorized operator. The factorization procedure essentially reduces to replacing the 
co-multiplicative operators comprising the initial factorization operator with some equivalent 
tridiagonal operators. The regularization parameter is determined by the maximum modulus of 
the eigenvalue of the metric tensor describing the curvilinear coordinate system. The non-linear 
terms are correctly taken into account in the scheme via an iterative process that is applied 
until a solution with a specified accuracy is obtained. At each iteration, a system of linear, 
algebraic equations with a tridiagonal matrix arises, which is solved numerically using a scalar 
fitting method. 

3 Computational results 

We have investigated the structure of the magnetic field in the accretion disk using the results 
of numerical simulations of MHD flows in the SS Cyg system [HJ |9]. The donor in this system 
is a red dwarf with mass 0.56 M Q , while the accretor is a white dwarf with mass 0.97 M Q . 
The orbital period of the system is P or b = 6.6 hr, and its semi-major axis is A = 2.05R & 
[46J. In our computations, the magnetic field at the surface of the white dwarf was specified 
to 10 5 G [4T| |4"8] . with the orientation of the magnetic axis given by 9 = 30° and <ft = 0°. 
The speed of gas flowing at the inner Lagrange point L\ was taken to be equal to the local 
sound speed, c s = 7.4 km/s, and the corresponding effective temperature of the donor to be 
4000 K. The gas density at L\ was taken to be p(Ll) = 1.1 x 10~ 7 g/cm 3 , and the mass-transfer 
rate to be M = 10~ 9 M /year. The solution was found in the region (— 0.56A < x < 0.56A, 
-0.56A <y< 0.56A, -0.28A < z < 0.28 A) on the geometrically adaptive grid [10] . 

The structure of the magnetic field in the accretion disk obtained in the numerical simula- 
tions is shown in Figs. [1] and |2] Figure [1] shows the distribution of the magnetic field B (in 
units of ^Airp(Li)Ail) in the equatorial (xy) plane of the disk for times of 12.75P or b (left) and 
13.36P or b (right). The light curve with arrows indicates the magnetic-field line passing through 
the point x = — 0A75A, y = 0.00 A The left-hand panel in Fig. [2] shows the radial distribution 
of the azimuth-averaged toroidal field in the disk B^ for the various times, while the right-hand 
panel in Fig. |2] shows the analogous distribution of the poloidal field B p = y ' B 2 + B 2 Z . 

Inspection of these figures indicates that three zones can be clearly distinguished in the 
disk: an inner zone where there is an intense generation of toroidal field due to the differential 
rotation of the disk, a zone of current sheets, and an outer zone of dissipation of the magnetic 
field. The radius of the inner zone is roughly 0.1A The current sheets are visible in Fig. [1] 
as dark rings. The magnetic- field lines change direction on either side of a current sheet. No 
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Figure 1: Magnetic-field distribution in the equatorial plane of the disk at times 12.75P or b (left) 
and 13.36P D rb (right). The light line with the arrows shows the magnetic field passing through 
the point x = -0.175 A, y = 0.00 A. 



fewer than three current sheets can be distinguished in the right-hand panel of Fig. [T] The 
current-sheet zone is located at distances from approximately 0.1^4 to 0.2 A. In Fig. [21 the 
current sheets correspond to points where the toroidal field B v vanishes. There is no intense 
generation of field by differential rotation in the outer zone, at distances exceeding 0.2v4, since 
the stellar magnetic field is very weak in this region. However, the right-hand panel of Fig. [2] 
shows that a weak poloidal field is generated in this zone. Overall, the magnetic field in the 
accretion disk is predominantly toroidal. 

The formation of current sheets occurs at the boundary of, or even inside, the inner zone. 
The change in the sign of B^ is associated with a variation in the rotation law in the disk. A 
transition zone where there is an intense loss of angular momentum forms near the star. In this 
zone, the angular velocity of the gas decreases from the Keplerian velocity to the rotational 
velocity of the stellar magnetic-field lines. This variation of the rotation law influences the 
generation of the toroidal field. 

These arguments can be clarified using the following picture. Consider the model depen- 
dence for the angular velocity u on radius r: 



u> 



u (r/r ) a 
l + (r/r ) a+3 / 2! 



(6) 



where r^ is some characteristic radius, approximately equal to the radius of the magnetosphere. 
We obtain when r -C r the asymptotic dependence u = u; (r/r ) a , which determines the 
rotation law in the region of the magnetosphere. In the opposite limit when r ^> r , expression 
(EJ) leads to a Keplerian rotation law, u = uo(r/r )~ 3 ^ 2 . Figure [3] illustrates the evolution of 
the magnetic field and formation of a current sheet. Figure [3] (first) shows the projection of the 
stellar magnetic-field lines onto the equatorial plane at the initial time. The subsequent panels 
show the pattern of the magnetic-field lines in the disk obtained by transforming the initial 
field pattern using the rotation law (jBJ) at times 27r/w (second), 47i/u (third) and 107i/u Q 
(fourth), respectively. Figure El shows that the current sheet and magnetic zones can already 
clearly form after the first several revolutions of the disk. However, the subsequent evolution of 
these structures can appreciably influence the dissipation of the magnetic field, radial motions, 
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Figure 2: 
magnetic 



Radial distribution of the azimuth-averaged toroidal B^ (left) and poloidal B p (right) 
fields in the disk for various times. 




Figure 3: Visual model for the formation of two magnetic zones and current sheets in an 
accretion disk. Schematics of the magnetic-field lines in the equatorial plane are shown for (left 
to right) the initial time and after one, two, and five revolutions of the disk. 

and dynamo processes. 

It follows that no current sheets should form in the case of equilibrium rotation of the 
star [3], when the corotation radius (distance where the rotational velocities of the field lines 
and of the disk material are equal) is equal to the magnetosphere radius. However, the disk 
material in such a star should be strongly wound up. Current sheets should form in systems 
where the corotation radius exceeds the magnetosphere radius. This condition is associated 
with the relationship between the accretor's rotational angular velocity, magnetic field, and 
accretion rate. Our computations assumed that the accretor rotated synchronously, so that 
this condition is obviously satisfied. Note that the generation of large-scale magnetic fields in 
galaxies by the dynamo mechanism together with the presence of a complex rotation law can 
also lead to the formation of magnetic zones separated by current sheets [27]. Such a current 
ring is present in the galaxy M31, for example, at a radius of approximately 3 kpc. 

The formed current sheet is carried to the outer part of the disk via decretion of the disk 
material. After some time, a new current sheet is formed in its place. Thus, several current 
sheets can exist simultaneously in the disk. Figure H] shows the distribution of the density 
(gray scale) and velocity (arrows) in the equatorial plane of the disk at times 12.75P or b (left) 
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Figure 4: Distribution of the density (gray scale) and velocity (arrows) in the equatorial plane 
of the disk at times 12.75P or b (left) and 13.36P or b (right). The lines with arrows show stream 
lines emerging from the points x = —0.175^4, y = 0.00A and x = —0.05A, y = 0.00A 



and 13.36P or b (right). The lines with arrows show stream lines emerging from the points 
x = — 0.05A, y = 0.00A (inner zone) and x = — 0.175A, y = 0.00A (outer zone). A comparison 
of the figure panels indicates that an alternation of accretion and decretion regimes occurs in 
the disk (especially its inner parts). The distribution of the azimuth-averaged radial velocity 
in the disk at various times is presented in Fig. El which clearly shows zones of accretion and 
decretion and their evolution with time. 

The alternation of accretion and decretion regimes in the inner parts of the disk can be 
understood if the magnetic field is quasi-periodic. In fact, the radial gradient in the pressure 
due to the toroidal magnetic field is increased by the generation of this field, ultimately causing 
the accretion to cease. This should lead to a flow of the field to outer parts of the disk, 
decreasing the magnetic pressure, and giving rise to a transition back to an accretion regime. 
The variations in the field distribution in the inner zone can be seen in Figs. [H [2j 

The generated magnetic field interacts intensely with spiral waves in the inner zone. Of 
most interest is the interaction of the magnetic field with a spiral precessional wave that can 
arise in the inner, unperturbed parts of the disk due to retrograde precession of the binary orbit 
[16]. In the gas-dynamical case, the increase in the radial flow of matter behind this wave leads 
to a growth in the accretion rate and the formation of a compact region of enhanced energy 
release at the surface of the accretor. In the magnetic case, the accretion acquires a column-like 
character. Consequently, the accretion rate should grow at times when the precessional wave 
approaches the stellar surface near the magnetic poles. In the coordinate system considered, 
the rotational period of the spiral wave is approximately equal to the orbital period of the 
system, if we disregard the longer precessional period. Thus, approximately twice per period, 
we should observe flares in the accretion rate, during the passage of the spiral wave by the 
North and South magnetic poles of the star. Figure [6] shows the evolution of the accretion rate 
onto the star over several orbital periods. The behavior shown is in good agreement with this 
picture. The maximum amplitude of fluctuations in the accretion rate is about 15%. Note that 
this amplitude can reach 200% in the case of ideal magnetic gas dynamics without diffusion 
effects [9]. The observed flux curves of SS Cyg in the UV between outbursts display similar 
behavior [H?] . with the flux fluctuations reaching 70%. 
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Figure 5: Radial distribution of the azimuth-averaged radial velocity in the disk v r for various 
times. 
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Figure 6: Evolution of the accretion rate onto the white dwarf. 



4 Generation of magnetic field 

4.1 Basic equations 

To elucidate the main regularities in the generation of magnetic field in the accretion disk, 
we consider a simpler model in which the magnetic-dipole axis of the star coincides with its 
rotational axis. In this case, we can assume axial symmetry when describing the magnetic field 
averaged over the turbulent pulsations and the azimuth [8] . The averaged magnetic field_| can 
conveniently be represented as B = B* + b, where B* is the magnetic field of the star and b 
the magnetic field generated by inductional currents in the accretion disk. The initial induction 
equation ([3]) taking into account the dynamo effect can be written 



db 

-=Vx(vxB, + vxb + ab- r?(V x b)) . 
at 



(7) 



1- We will denote the averaged magnetic field in this section B. This should not lead to confusion, since the 
initial, unaveraged field that satisfies ([3]) , is not used anywhere in this section. 



Since the magnetic field is poloidal, (B^ = 0, B* p = B*), the toroidal component of the 
generated field is b^, = B^,. We denote b p = V x (An v ), where n^, is a unit vector in the 
azimuthal direction. Using this notation, Eq. ([7]) can be rewritten as the system 

8 A 1 / A \ 

— + -(v p ■ V)(M) = aB v + V (v 2 A - -J , (8) 

^ = r(B p ■ VV - rV • (^v p ) - \va • V(M) - a ( V M - ^) + r, (v 2 B, - ^) , (9) 

where u = v v /r is the rotational angular velocity in the disk. 

Under the conditions in the accretion disk, these equations can be appreciably simplified. 
Since |u r |, \v z \ <C \v v \ in the accretion disk, we can set v p = 0. Moreover, we can omit the 
negligibly small third and fourth terms in the right-hand side of (j^J), which correspond to the a 2 
dynamo. The radial derivatives can be neglected in the diffusion terms compared to the vertical 
derivatives, since their ratio is z/r <C 1. Further, using the estimates dco/dz ~ z/rdco/dr, 
\b z du)/dz\ <^C \b r duj/dr\ and 

B„ = -B a l-±) -, B* z = --B a l-^J , (10) 

we find 

z dA 

r{B p -V)u = -gB -g—, (11) 

where B a is the magnetic field at the stellar surface, g = rdu/dr a measure of the differential 
rotation, and B = B a (R a /r) 3 the characteristic stellar magnetic field in the disk. Finally, we 
will assume that g and rj depend only on r. 
These simplifications lead to the system 



dA d 2 A 

-^^otB. + n- 



aB v + r]—, (12) 



9B V z dA d 2 B^ 

^r = -r 9B °- 9 ^ +r] ^^- (13) 

These equations have a simple physical meaning. Equation ( TT2|) describes the evolution of 
the poloidal magnetic field in the disk. The first term on the right-hand side determines the 
generation of poloidal field from toroidal field via the dynamo effect, while the second term 
determines the dissipation of the poloidal magnetic field. Equation (1T3|) describes the evolution 
of the toroidal magnetic field in the disk. The first and second terms on the right-hand side 
of this equation determine the generation of toroidal field from the stellar magnetic field and 
the generated poloidal field due to the differential rotation of the disk. The last term on the 
right-hand side describes the dissipation of the toroidal magnetic field. 

4.2 Generation of field by differential rotation 

The dynamo effect is weak in the inner parts of the disk, and magnetic field is generated by 
differential rotation in the disk. In this case, we can neglect the aB^ term in (IT2"|) . As a result, 
we find A = 0, so that only toroidal magnetic field will be generated. The remaining Eq. ( 1131) 
acquires the form 

dB v z d 2 B v 

-W = r 9B ° + ^ (14) 
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Figure 7: Vertical and radial profiles of the toroidal magnetic field in the inner parts of the 
accretion disk. 



The general solution of (JHJ), satisfying the initial condition B v {t 
0) = B v (z = H) = can be written [50] 



0) = and the boundary 



conditions B v (z 



B u 



25 Vi0 Yl 



n=l 



:-iy 



[Tvn] 



-(im) 2 t/t d 



f-Knz 

sm hr 



(15) 



where B^q = gB Htd/r, and td = H 2 /rj is the diffusion time. In early stages of the field 



generation when t -C td, the toroidal field grows linearly, B v - 
describes the generation of field in the ideal MHD case with 77 



gB zt/r. This solution also 
= 0. Note, however, that this 



leads to a discontinuity in B v at the surface of the disk at z = H . This means that, near this 
surface, an important role is played by dissipation of the field, so that this cannot be neglected, 
at least in this region. In the other limiting case when i > irf we obtain the stationary solution 



B u 



6 Bv '°H 



1 



z 

1p 



(16) 



The maximum magnetic field B v is reached at the point z = H/\/3. The mean field averaged 
over height in the disk is B v = B^/24. 

The coefficient of B^q in (TIB"]) describes the radial structure of the generated toroidal field. 
Let us consider the model dependence ([6]) of the rotational angular velocity in the disk. If we 



„a-2 



assume that H oc r and r\ oc r, we can obtain the asymptotic relations B^p oc r a z as r < r 



and B u 



-7/2 



'<p,o oc r" v " as r ^> r . The behavior B^q as r — > depends appreciably on the parameter 
a, which characterizes the rotation of the disk in the region of the stellar magnetosphere: 
B^q — > 00 as a < 2, B^q — > const as a = 2, and B^q — > as a > 2. 

The left-hand panel in Fig. [7] shows characteristic vertical profiles of the stationary toroidal 
magnetic field in the accretion disk for a = 3 at various distances from the star. The right- 
hand panel of Fig. [7] shows the characteristic radial profiles of the field for various values of 
a. The sign of the magnetic field changes at the point r = (2a/3) 1/, ^ a+3//2 Vo. This means that 
neighboring magnetic-field lines will have opposite directions at this distance, leading to the 
formation of a current sheet separating two magnetic rings. The magnetic fields in the inner 
and outer magnetic rings have opposite signs. Recall that the numerical simulations lead to a 
similar picture (Fig. [2]). 
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Figure 8: Illustration of the mechanism for the generation of a quadrupolar component of the 
toroidal field in the accretion disk in the case of an inclinedmagnetic moment. 

The antisymmetric character of the dependence of B v on z is related to the assumption that 
the magnetic and rotational axes of the star coincide. The antisymmetry of the field about the 
equatorial plane is disrupted if the magnetic axis is inclined. This situation is illustrated in 
Fig. [SI a clear vertical gradient of the field (shown by the light arrows) arises in the disk when 
the magnetic axis is inclined. In the right-hand panel, the field increases from top to bottom, 
and in the left-hand panel from bottom to top. As a result, primarily the stronger field will be 
amplified. Taking into account the rotation direction in the disk, the field directed toward the 
observer will grow more rapidly in the right-hand panel, and the field directed away from the 
observer in the left-hand panel. In both cases, the resulting toroidal field is positive. This type 
of toroidal-field distribution is also obtained in the computations. This effect will be expressed 
most strongly when the inclination of the magnetic axis is 9 — 45°. This reasoning bring us 
to the following conclusion. An antisymmetric toroidal field (dipolar mode) will be generated 
only when 9 = 0° and 90°. Any small deviation from these values will cause the quadrupolar 
mode to suppress the dipolar mode. Therefore, the generation of a purely dipolar mode in the 
accretion disks of intermediate polars may be unstable. The quadrupolar mode is more stable. 
Recall also that the dynamo generation of magnetic field in the outer part of the disk also leads 
primarily to the growth of the quadrupolar field component |27| . 
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The dynamo effect can play an important role in the generation of magnetic field in the outer 
parts of the disk. Therefore, we must use the full system of equations (fi~2| [T3|) to describe the 
magnetic field. The quantity a, in these equations is determined by the mean helicity of the 
velocity fluctuations 5v. 

(17) 
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Figure 9: Distribution of the mean helicity a in the accretion disk at various distances from 
the accretor. 



where r is the characteristic correlation time. The angular brackets denote an average over 
the ensemble of turbulent pulsations and over azimuth. In both averages, the characteristic 
correlation time can be estimated as r = H/(aTC s ). 

Figure [9] shows vertical distributions of a for non-azimuthal velocity fluctuations in the 
accretion disk at various distances from the accretor obtained via three-dimensional MHD 
numerical simulations. The helicity possesses a mirror antisymmetry about the plane of the 
disk, which is a necessary condition for the operation of the dynamo mechanism [27] . Similar 
profiles are probably obtained for turbulent velocity fluctuations as well. To simplify our further 
analysis, we will use values of a averaged over half the disk height. 

We will restrict our analysis to stationary solutions of (fT2" |[T5|) . In this case, these equations 
can be reduced to a single equation for the toroidal field: 

d 3 B u 



dz 3 



ag 
rf 



+ -jB v 



gBo 

7/r 



0. 



\i 



-> -B u 



This equation is invariant to the transformation z — > —z simultaneous with B^ 
This means that it describes the generation of dipolar (antisymmetric) magnetic field, for 
which B v (—z) = —B v (z). As was shown in the previous section, this type of field is excited 
by the stellar field subject to the differential rotation of the disk. However, the quadrupolar 



(symmetric) field, for which B v {—z) 
the star, where gB H 3 //(r]r) < 1 [271 



B v {z) 



may become dominant at large distances from 



12 



Due to the linearity of f JT8|) . a general solution can be written in the form B, 
where B^ 
equations 



M 



>(?) 



R K ' J- R y 

and BJp' describe the dipolar and quadrupolar field components and satisfy the 
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These equations can be rewritten in dimensionless form. We denote B 
where ( = z/H. Equations ( FT9| 120]) then take the form 

//' + Dfa + 1 = 0, 



(«»,«) 



fa" + £>/, 



0. 



(19) 
(20) 

B<f,flfd,q{0) 

(21) 
(22) 



where D = agH 3 /rj 2 is the dynamo number and a prime denotes differentiation with respect to 
(. These equations can be solved using the boundary conditions fa(0) = /d(l) = 0, /</'(l) = — 1) 
/ g '(0) = f q (l) = fq'iX) = 0. The dynamo number D should be negative in outer parts of the 
disk when z > since the differential-rotation measure g is negative in these regions. 
The solution of ([21]) for the dipolar mode can be written 



f d = ^ + C d %^ + e~^ 



Cf sin 



'v^iV 



C ) + Cf cos 



'V3N 



(23) 



where iV = IDI 1 ' 3 . The constants of integration C d ' ' in this expression can be found from 
the boundary conditions. In the limit of small dynamo numbers as iV — > 0, we obtain the 
asymptotic relation fa = £(1 — ( 2 )/6 + O(N), corresponding to the solution ffT6|) . 
The solution of (|2"2"1) for the quadrupolar mode can be written [2] 



fq — C„ 



MC-D _ 2e -T(C-D 



cos 



'v^A^ 



(C 




(24) 



where C q is the amplitude. The quantity iV satisfies the relation 

, V3N\ _3N 

2 cos I — - — I + e 2 



(25) 



The roots of this equation determine the discrete set of allowed values of N. The first three 
roots are JVi = 1.85, N 2 = 5.44, iV 3 = 9.07. 

The obtained dependences ( 123]) and (1241) determine the vertical structure of the field. Note 
that, when the dynamo effect is included, the radial structure of the field is described by -8^,0, 
as before. The quadrupolar component will probably always dominate in the inner magnetic 
ring when the magnetic axis is inclined. The quadrupolar field can also become dominant in 
the outer magnetic ring with time, since this component grows more rapidly than the dipolar 
field due to the dynamo effect [27]. This magnetic-field structure in the inner and outer rings 
is clearly demonstrated by our computations. 
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Figure 10: Distribution of the density and velocity in the xz plane. 

5 Conclusion 

We have investigated the magnetic-field structure in accretion disks in semi-detached binary 
systems based on three-dimensional numerical simulations. These assumed that the intrinsic 
magnetic field of the accretor is dipolar, with the dipole axis, in general, inclined to the rotation 
axis. The numerical simulations took into account the diffusion of magnetic field due to the 
dissipation of currents in turbulent vortices and magnetic buoyancy. We considered the SS Cyg 
system as an example. 

The results of the simulations show that the disk magnetic field is predominantly toroidal. 
Three zones can be distinguished in the disk: an inner zone with the intense generation of 
toroidal field via the differential rotation of the disk, a zone with current sheets, and an outer 
zone of dissipation of the magnetic field. A current sheet forms in the inner zone, due to the 
variation of the rotation law for the disk material near the stellar magnetosphere. Further, 
this current sheet is carried into the outer part of the disk with decreted material, with a new 
current sheet forming in its place after some time. Therefore, several current sheets can be 
present simultaneously in the intermediate zone. The alternation of accretion and decretion 
regimes in the inner parts of the disk is probably due to the generation of magnetic field and 
the associated increase in the corresponding magnetic-pressure gradient, which acts to stop the 
accretion. This is followed by an outflow of the field into the outer part of the disk, decreasing 
the magnetic pressure and bringing about a transition back to an accretion regime. 

Current sheets should form in systems in which the corotation radius exceeds the mag- 
netosphere radius. This condition is determined by the relationship between the accretor's 
rotational angular velocity, magnetic field, and accretion rate. The field dissipation and tem- 
perature are enhanced in current sheets, leading to the emission of UV and X-ray radiation. 
Moreover, the total pressure is reduced in current sheets due to the decrease in the magnetic 
field. Therefore, the disk thickness should decrease in the vicinity of current sheets, as is 
confirmed by the numerical simulations. Figure [TD] shows the density distribution in the xz 



14 



plane in the disk at time 13.36P or b. The disk thickness decreases slightly in the intermediate 
zone (0.1 < r/A < 0.2), where the current sheets are located (Fig. [T]), compared to the inner 
and outer zones. Note that, in the purely gas-dynamic case [5] the disk thickness increases 
monotonically with distance. 

The magnetic field in the inner zone interacts intensely with a spiral precessional wave. As 
a result, the accretion rate grows when the precessional wave approaches the stellar surface 
near the magnetic poles. This gives rise to a flare in the accretion rate (with an amplitude of 
the order of 15%) twice per period, when the spiral wave passes the North and South magnetic 
poles of the star. This behavior of the accretion rate is in good agreement with the observed 
UV flux variations of SS Cyg between outbursts. 

A vertical gradient of the field arises in the disk when the magnetic axis is inclined. There- 
fore, the generation of toroidal field in the disk acquires a symmetrical (quadrupolar) character 
about the equatorial plane. This effect is most strongly expressed when the inclination of the 
magnetic axis is 9 = 45°. An antisymmetric toroidal field (dipolar mode) will be generated 
only when 9 = 0° and 90°. Small deviations from these values lead to a suppression of the 
dipolar mode by the quadrupolar mode. Thus, the generation of a purely dipolar mode in the 
accretion disks of intermediate polars may be unstable. The quadrupolar mode is more stable. 

Our computations show that motion in the accretion disk displays a mirror antisymmetrical 
helicity. Therefore, the magnetic field in the outer parts of the disk may be partially generated 
by the dynamo mechanism. Here, the quadrupolar mode again grows most rapidly. However, 
the role of the dynamo is small, since dynamical effects dominate in the disk. 
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